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SUMMARY 

A  theoretical  as  well  as  an  experimental  investigation  of  a 
plane,  two-dimensional,  turbulent  wall  jet  and  a  turbulent  wall  jet  be¬ 
neath  a  secondary  uniform  stream  is  presented. 

For  the  case  of  the  plane  turbulent  wall  jet,  the  growth  of  the 
jet  as  well  as  the  maximum  velocity  decay  are  predicted  based  on  the  as¬ 
sumption  that  the  shear  stress  distribution  as  well  as  the  velocity  dis¬ 
tribution  across  the  jet  remain  simi  lar  over  the  useful  range  of  downstream 
positions.  A  suggested  form  for  the  velocity  distribution  is  compared  with 
the  experimental  results  of  this  investigation  as  well  as  with  the  results 
of  several  previous  investigators. 

The  approach  used  in  the  analysis  of  the  plane  turbulent  wall 
jet  is  then  extended  to  the  more  general  case  of  a  turbulent  wal  I  jet 
beneath  a  secondary  uniform  stream.  The  velocity  distribution  and  shear 
stress  distribution  are  expressed  in  terms  of  a  number  of  similarity  func¬ 
tions  and  the  ratio  of  the  free  stream  velocity  to  the  maximum  velocity 
in  the  jet  profi le.  Forms  for  these  simi larity  functions  are  also  sug¬ 
gested.  The  growth  of  the  jet  and  decay  of  maximum  velocity  are  again 
predicted  based  on  the  results  obtained  from  the  plane  turbulent  wall  jet 
analysis.  The  experimental  results  of  this  investigation  and  the  results 
of  previous  investigators  are  compared  with  the  velocity  profiles  as  pre¬ 
dicted  by  this  analysis. 
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NOMENCLATURE 

A 

constant  equal  to  1.4794 

B 

constant  equal  to  .67753 

C 

function  of  jet  Reynolds  number 

D 

function  of  jet  outlet  conditions  and 

jet  Reynolds  number 

F 

function  of  t 

M 

function  of  jet  outlet  conditions  and 

jet  Reynolds  number 

S 

shear  similarity  function  of  1] 

T 

shear  similarity  function  of  "f) 

U 

mean  velocity  in  x  direction 

V 

mean  velocity  in  y  direction 

w 

shear  simi  larity  function  of  T| 

*  1 

shear  simi  larity  function  of  T) 

d 

jet  slot  height 

e 

exponential  function 

f 

simi  larity  function  of  T1 

g 

simi  larity  function  of  T) 

h 

simi  larity  function  of  T] 

j 

simi  larity  function  of  T1 

u 

non-dimensional  x  velocity  component. 

U/Vj 

V 

non-dimensional  y  velocity  component. 

V/Vj 

X 

longitudinal  distance  from  jet  exit 

•  y 

t 

vertical  distance  above  wall 

a  constant,  fa(ipdT)  (also  constant  in  Glauert's  analysis) 

P  constant,  P  f  (T)> C2g <T|)  -  I  ]dTl 

£>  CD 

V  constant,  P  g(T))[g(T))  -  IldT)  +  2  J  f  (Tph(T]) d*[) 

o 

6  characteristic  jet  thickness 

Q  velocity  ratio,  U^/Vj 

T]  similarity  variable,  y/6 

\  simi  larity  parameter 

p,  viscosity  coefficient 

u  kinematic  viscosity 

p  mass  density 

00  CO 

a  constant,  J  h(Tp[2g(Tl)  -  I  ]dT)  +  2  J  f  (T|)  j (Tl) dTl 

o  o 

t  turbulent  shearing  stress 
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U  shearing  velocity,  X 

*  V  p 

u*  non-dimensional  shearing  velocity,  U*/Vj 

Vj  jet  outlet  velocity 

x  non-dimensional  longitudinal  coordinate,  x/d 

y  non-dimensional  vertical  coordinate,  y/d 

Umax  maximum  velocity  in  jet  profile 

Um  non-dimensional  maximum  velocity,  Umax^V j 
6  non-dimensional  jet  thickness,  6/d 

t w  wall  shearing  stress 

UT  skin  friction  velocity, 

uT  non-dimensional  skin  friction  velocity,  UT/Vj 


I  V 


Cf  skin  friction  coefficient 

T)#  value  of  T|  where  S(T|)  vanishes 

T]m  value  of  1)  where  f*(T|)  vanishes 
Re0  jet  outlet  Reynolds  number 

Kj  constant  (i  =  I,  2,  •••  ) 

4»  free  stream  velocity 

T|i  /s  function  of  C/Um 


I  INTRODUCTION 


The  behavior  of  a  jet  blowing  tangentially  to  a  solid  surface 
is  of  great  i nterest  because  of  the  application  of  such  jets  to  boundary 
layer  control  of  aircraft.  In  the  absence  of  a  pressure  gradient,  the 
flow  field  produced  by  such  a  jet,  commonly  called  a  wall  jet,  is  a  com¬ 
bination  of  a  boundary  layer  and  a  free  turbulent  flow.  This  flow  field 
may  be  further  complicated  by  the  presence  of  a  secondary  free  stream 
as  in  the  case  of  any  normal  boundary  layer  control  application. 

For  the  case  of  a  plane  wall  jet,  Glauert  (I)  has  obtained  a 
theoretical  solution  of  the  boundary  layer  equations  for  both  laminar  and 
turbulent  flow.  However,  the  analysis  of  the  laminar  wall  jet  in  which 
a  similarity  solution  was  found  is  of  little  practical  interest  because 
of  the  extremely  low  critical  Reynolds  number  of  mixing  layer  flows. 

For  the  turbulent  case,  a  near  simi larity  solution  was  obtained  by  as¬ 
suming  that  the  variation  of  shear  stress  with  mean  velocity  near  the 
wall  was  the  same  as  that  in  a  turbulent  pipe  flow.  Outside  of  the  ve¬ 
locity  maximum,  an  eddy  viscosity  which  remained  constant  across  the 
breadth  of  the  jet  was  assumed,  the  value  being  chosen  to  give  the  best 
agreement  with  experiment.  As  a  result  of  these  assumptions,  the  shear 
stress  at  the  maximum  velocity  location  necessari ly  vanished,  a  result 
which  has  been  questioned  by  several  recent  experimental  investigators 
(2,  3,  4). 

The  more  practical  case  of  a  wall  jet  beneath  a  uniform  free 
stream  is  not  as  easi ly  treated  by  the  concept  of  simi lar  profi les  since 
the  ratio  of  tree  stream  velocity  to  maximum  velocity  in  the  jet  profi le 
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varies  continuously  with  the  downstream  coordinate.  However,  tentative 
theories  in  which  a  form  of  downstream  similarity  is  assumed  have  been 
suggested  by  George  (5)  and  Carriere,  Ei che Ibrenner,  and  Poi sson-Qui nton 
(6). 

As  in  the  case  of  the  plane  turbulent  wall  jet,  the  correspon¬ 
dence  of  the  point  of  zero  shear  stress  with  the  maximum  velocity  loca¬ 
tion  has  also  been  questioned  experimentally.  This  apparent  lack  of 
dependence  of  the  Reynolds  stress  on  the  mean  velocity  gradient  seems  to 
disprove  the  validity  of  mixing  length  theory  and,  in  particular,  the 
concept  of  an  eddy  viscosity  when  applied  to  an  asymmetric  velocity  distri¬ 
bution.  Furthermore,  Glauert's  assumption  that  the  flow  may  be  divided, 
for  purposes  of  calculation,  into  layers  each  having  the  properties  of 
some  known  type  of  shear  flow  is  also  subject  to  question. 

Based  upon  the  experimental  results  of  this  investigation  as 
well  as  those  of  George  (7)  and  Myers,  Schauer,  and  Eusti s  (8)  for  the 
plane  turbulent  wall  jet,  the  present  analysis  assumes  that  similarity 
exists  with  respect  to  the  velocity  profiles  over  the  entire  useful  range 
of  downstream  positions.  The  assumption  that  the  shear  stress  distribu¬ 
tion  also  remains  similar  leads  to  a  relation  between  the  jet  thickness 
and  maximum  velocity  in  the  jet  profile.  This  relation,  together  with 
the  consideration  of  the  effect  of  introducing  a  "wall"  into  a  two- 
dimensional  turbulent  jet,  yields  an  approximate  relation  for  the  growth 
of  the  jet  and  decay  of  maximum  velocity  in  the  plane  turbulent  wall  jet. 

The  present  analysis  of  the  turbulent  wall  jet  beneath  an  ex¬ 
ternal  free  stream  expresses  the  velocity  profi les  and  also  the  shear 
stress  distribution  in  terms  of  a  series  of  similarity  functions.  The 


results  for  the  plane  turbulent  wall  jet  are  appropriately  modified  in 
order  to  predict  the  growth  of  the  jet  and  decay  of  maximum  velocity  in 
the  presence  of  a  free  stream.  Because  of  the  difficulties  encountered 
in  obtaining  precise  analytic  representations  for  the  similarity  functions, 
an  approximate  form  for  the  velocity  distribution  which  includes  only  four 
such  functions  has  been  adopted.  This  approximation  remains  valid  as  long 
as  the  maximum  velocity  in  the  profi le  i s  at  least  as  great  as  the  free 
stream  velocity. 
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II  THEORETICAL  ANALYSIS 

Part  I :  Plana  Turbulent  Wal I  Jet 

A.  Equations  of  Motion 

The  turbulent  wall  jet  may  be  regarded  as  a  combination  of  a 
boundary  layer  and  a  free  turbulent  flow.  Since  problems  in  free  turbu¬ 
lent  flow  are  of  a  "boundary  layer"  nature  (i.e.,  the  transverse  gradients 
are  large  and  the  region  in  which  a  solution  is  sought  does  not  extend 
far  in  the  transverse  direction  as  compared  with  the  longitudinal  direc¬ 
tion  cf  flow),  it  is  expected  that  the  conventional  boundary  layer  equa¬ 
tions  are  also  applicable  to  the  case  of  b  wall  jet.  For  a  two-dimensional, 
incompressible  flow  in  the  absence  of  a  pressure  gradient,  these  equations 
are  (9) 

U  ail  +  V  I  =  i  ii  Cl.  la) 

dx  dy  p  dy 

and 

2U  +  SI  =  0  (l.lb) 

dx  dy 

where  U  and  V  refer  to  mean  velocities  and  t  denotes  the  turbulent  shear¬ 
ing  stress. 

Defining  the  shearing  velocity  as 

(1.2) 

(1.3a) 

(1.3b) 


U*;/T 


and  non-dimensional izing 


u  =  i l 

Vi 

V  =  ! 

Vj 

u#  = 

VJ 


( I .3c) 
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x  =  x  (1.3d) 

d 

y  =  £  (I .3e) 

d 

where  Vj  is  the  jet  outlet  velocity  and  d  the  jet  slot  height,  the  boun 

dary  layer  equations  become 

u  iu  +  v  Ju  =  (1.4a) 

3x  9y  3y 

and 

au,  +  av  =  o  ( i .  4b) 

3x  Sy 


B.  Simi lar  Velocity  Profi les 

For  the  plane  turbulent  wall  jet,  the  general  form  of  the  mean 
velocity  profile  is  as  shown  below,  the  maximum  velocity  in  the  profile 
being  denoted  by  Umax.  In  general,  the  flow  in  the  outer  part  of  the 

profile  is  virtually  the 
same  as  that  i n  a  two- 
dimensional  free  turbulent 
jet  provided  the  jet  mo¬ 
mentum  i s  very  much  greater 
than  the  wall  friction. 
Except  for  the  loss  of  mo- 
montum  due  to  this  surface 
friction,  the  flow  is  at 
constant  momentum  so  that 
at  least  the  overall  chai — 
acteri sties  of  the  flow  will 
depend  on  this  momentum. 
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|f  these  profiles  are  to  be  independent  of  the  downstream  co¬ 
ordinate  (i.e.,  if  there  i s  to  be  a  similarity  solution),  then  the  gen¬ 
eral  form  of  the  velocity  distribution  must  be 

U  =  f  (T))  (1.5) 

Umax 

where 


T\  =1  , 
6 


(1.6) 


6  being  some  characteristic  thickness  of  the  jet  profile.  Several  pro¬ 
perties  of  this  similarity  function  f(T])  can  be  readily  deduced,  namely 

f (0)  =  0  (1.7a) 

f  (eo)  =  0  (1.7b) 

and 


occurring  at  some  T]  =  T^. 

Assuming  that 


f 


max  = 


6  a  £ 
d 


(1.7c) 


(1.8) 


and 

Um  B  Umax  (1.9) 

Vj 

are  functions  only  of  x,  then  integration  of  the  continuity  equation  re¬ 
sults  in 

v  =  Um  dL  71  f(Tl)  -  A.  (5  Um)  J11  f(z)dz  (I.  10) 

dx  dx  o 

With  regard  to  the  characteristic  jet  thickness  6,  two  possible 
choices  are  available,  one  being  the  value  of  y  where  the  maximum  velocity 
occurs  and  the  other  being  that  value  of  y  in  the  outer  portion  of  the 


jet  where  the  velocity  Is  half  the  maximum  value.  Since  the  latter  choice 
is  more  easily  measured  experimentally,  it  is  usually  taken  as  the  repre¬ 
sentative  thickness  of  the  wall  jet.  This  then  imposes  another  condition 
on  the  similarity  function,  namely, 

f  C I )  =  1/2  ( 1 .  1 1 ) 


C.  Shear  Simi  laritv 

Making  use  of  equations  (1.5)  and  (1.10),  the  momentum  equation 
( I . la)  may  be  written 

(1.12) 


Um  dum  fa<Tp  -  Uffl  (6Um)  f'(Tp  J  f (z) dz  = 
dj?  6  dx  °  dy 


au*2 


Integrating  this  equation  with  respect  to  y  from  zero  to  °°  gives 

00 

_1_  la  -  -  Um  [26  dUtri  +  um  4  ]  J  f2  (11)  dT)  ( 1 . 13) 

w  .a  -  dx  n 


V.i-  P 


dx 


where  tw  is  the  wall  shearing  stress.  Defining  the  skin  friction  velo¬ 
city  as 


with 


and  letting 


”  0 


o  =  J  f2(Tl)dT) 


(1.14) 


(1.15) 


(1.16) 


equation  (1.13)  may  be  written 


3  _ 


a  A.  (6  Um2) 
dx 


Defining  a  skin  friction  coefficient  by 


:f  = _ i_-  2H 

u 


3 

T?W 

tn 


(1.17) 


(1.18) 
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then  from  equation  (1.17) 

Cf  *  -  20.  4,  <6  U  8)  (1.19) 

Um3  dx 

On  the  basis  of  Forthmann's  (10)  calculations,  it  may  be  as¬ 
sumed  that  si  mi larity  also  exists  with  respect  to  the  shear  stress  dis¬ 
tribution  in  a  turbulent  wall  jet  flow.  That  is, 

u*2  =  uT2  S(Tp  (1.20) 

where  S (T|)  is  the  shear  similarity  function.  This  is  also  consistent 
with  the  claim  of  Eskinazi  and  Kruka  (4)  that  it  is  a  characteristic  shear¬ 
ing  velocity  and  not  a  characteristic  mean  velocit'  which  brings  similar¬ 
ity  to  turbulent  velocities  (i.e.,  to  turbulent  shear  stress). 

Making  use  of  equation  (1.20),  the  momentum  equation  (1.12) 
may  then  be  written 

m  3 

um  f S op  -  iJn  Jl_  (6  Um)  f(Tl)  J  f(z)dz  «  S'(T1) 

dx  6  dx  °  5 

The  conditions  which  must  be  satisfied  by  the  similarity  function  S(T)) 
are 


S(0)  =  I 
S(»)  =  0 
S'(0)  =  0 

Making  use  of  equation  (1.17) 


( 1 .2  la) 
( 1 .2 1 b) 
(l.2lc) 

the  momentum  equation  may  be 


written 

Um  [  f 2 (71)  -  f '  (T|)  J  f(z)dz  +  2a  S'(T))  ]  = 
dx  ° 

T], 

Hm_d£[  f'ttj)  J  f (z)dz  -  as '(7))  ]  (1.22) 

5  dx  o 

If  there  i s  to  be  a  si  mi larity  solution  (i.e.,  a  solution  which  is  inde¬ 
pendent  of  x),  then  both  sides  of  this  equation  must  vary  in  the  same 
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manner  with  x.  This  implies,  for  some  constant  \, 

Hal  di  =  -  X  Um  d^jm 
5  dx  dx 

or 

dUm  +  i  =0  (1.23) 

Um  dx  5  dx 

Therefore,  the  momentum  equation  becomes 

jn 

f3  (Tl)  f  (X  -  I)  f'(Tl)  J  f  (z)  dz  +  (2  -  X)aS '  (T|)  =  0  (1.24) 

o 

which  is  the  similarity  equation  governing  the  flow  in  a  plane  turbulent 
wa  I  I  jet. 

The  general  shape  of  the  shear  function  is  shown  below 


Denoting  the  value  of  T]  where  S(T))  vanishes  by  T|  ,  then  equation  (1.24) 
may  be  integrated  giving 


10- 


(x  -  i)  f(V 


f  C ”H )  dTJ  +'  (X  -  2) 


J*  fs  (7))dT)  =  0 
% 


(1.25) 


D.  Simi laritv  Function 

A  solution  of  equation  (1.24)  cannot  be  found  unless  some  fur¬ 
ther  assumptions  are  made  concerning  the  velocity  distribution  f (T|)  or 
the  shear  distribution  SO)).  Glauert  (I)  has  obtained  a  solution  for  the 
plane  turbulent  wall  jet  assuming  that  the  inner  part  of  the  flow  near 
the  wall  is  governed  by  Bias! us'  formula  and  the  outer  part  by  Prandtl's 
hypothesis,  the  matching  taking  place  at  the  maximum  velocity  point. 
Blasius'  formula  implies  that  U  ~  y1^7  near  the  wall  while  Prandtl's  hy¬ 
pothesis  states  that  the  eddy  viscosity  is  constant  across  the  outer  por¬ 
tion  of  the  flow  and  is  proportional  to  the  product  of  the  maximum  mean 
velocity  and  the  thickness  of  the  jet. 

As  a  result  of  the  assumptions  concerning  the  eddy  viscosity, 
the  solution  involves  one  disposable  constant,  and  exact  similarity  of 
the  velocity  profiles  at  all  downstream  positions  does  not  exist  since 
the  eddy  viscosity  in  the  inner  and  outer  parts  of  the  jet  vary  in  a 
slightly  .different  manner  with  the  Reynolds  number  based  on  jet  thickness 
&nd  maximum  velocity.  However,  the  velocity  profiles  are  relatively  in¬ 
variant  for  given  orders  of  the  Reynolds  number,  the  general  nature  of 
these  profiles  being  characterized  by  a  parameter  o  which  Is  indicative 
of  the  proportional  extent  of  the  inner  layer,  between  the  wall  and  point 
of  maximum  velocity,  to  that  of  the  outer  layer,  between  the  velocity 
maximum  and  the  stationary  fluid  far  from  the  wall. 

The  relative  invariance  of  the  velocity  profiles  in  a  plane 
turbulent  wall  jet  has  been  adequately  verified  by  previous  experimental 


investigations  as  well  as  by  the  present  one.  Those  investigators  who 
have  related  their  results  to  Glauert's  solution  have  all  adopted  pro¬ 
files  corresponding  to  a  narrow  range  of  the  parameter  a  (i.e.,  1.2  < 
a<  1.3).  This  possibly  indicates  that  complete  similarity  does  exist 
to  a  greater  extent  than  predicted  by  Glauert's  assumptions  concerning 
the  eddy  viscosity  which  are  of  doubtful  validity. 

On  the  basis  of  this  information,  it  will  therefore  be  assumed 
that  the  velocity  profiles  in  a  plane  turbulent  wall  jet  can  be  adequate¬ 
ly  represented  by  a  universal  similarity  function  f(Tp.  The  shear  stress 
distribution  is  then  given  by  the  simi  larity  function  S (T))  which  is  ex- 
pressable  in  terms  of  f (T))  by  means  of  equation  (1.24). 

From  the  standpoint  of  obtaining  some  sort  of  solution  of  equa¬ 
tion  (1.24),  it  is  much  more  practical,  although  not  as  mathematically 
satisfying,  to  make  some  plausible  assumptions  concerning  the  velocity 
distribution,  since  the  general  shape  of  the  mean  velocity  profile  in  a 
plane  turbulent  wall  jet  is  reasonably  well  defined.  One  such  assumption 
is  that  the  inner  portion  of  the  flow,  which  closely  approximates  a  tur¬ 
bulent  boundary  layer,  is  related  to  Blasius'  formula  which  implies  f (T|)  ~ 
T]1^7  near  the  wall. 

For  the  case  of  a  free  jet  boundary,  H.  Goertler  has  obtained 
a  solution  for  the  velocity  distribution  in  terms  of  the  Error  function 
making  use  of  Prandtl's  hypothesis.  Therefore,  a  reasonable  form  for 
the  similarity  function  f(T|)  suggested  by  Blasius'  formula  and  Goertler's 
solution  is 

f(Tl)  =ATll/7  [  I  -  erf  (BT|)  ]  (1.26) 
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where  A  and  B  are  constants.  This  function  automatically  satisfies  con¬ 
ditions  (1.7a)  and  (1.7b)  so  that  the  remaining  conditions  (1.7c)  and 
(1.11)  serve  to  determine  the  constants  A  and  3.  These  constants  were 
found  to  be  1.4794  and  .67753  respectively  so  that 

f(7])  =  1.4794  C  I  -  erf  (.67753  7])  ]  (1.27) 


with  71  =  .  16539. 


From  equation  (1.24) 

T\ 

S' (7))  =  f3 (Tl)  +  X  -  I  f  •  (71)  f  f  (z)dz 

a(X  -  2)  a(X  -  2)  o 


( I .28) 


which  can  be  integrated  by  parts  to  give 

„7]  7) 

S (71)  *  I  +  X  -  I  f (71)  f  f (z)dz  “if  f3(z)dz  (1.29) 
a(X  -  2)  t  a  * 

since  S(0)  =  I.  One  obvious  shortcoming  of  the  assumed  form  of  f(7))  is 
the  fact  that  f '  (71)  increases  without  limit  as  7)  approaches  zero.  How¬ 
ever,  since 


Tl 

I  im  f '  (71)  f  f  (z)dz  =  0 
THO 


S'(71)  as  defined  by  equation  (1.28)  exhibits  no  such  behavior. 

The  function  f(Tp  was  programed  for  the  computer  and  the  values 

71  71 

of  f'(Tl),  f(Tp,  £  f(z)dz,  and  J  f3(z)dz  are  tabulated  in  Table  I.  The 
procedure  used  in  dealing  with  the  singularity  f * (0)  is  described  in  the 
Appendix.  The  constant  a  defined  by  equation  (1.16)  was  found  to  be 


00 

a  =  f3  (TpdTl  =  .75437 


(1.30) 


as 


From  equation  (1.25),  the  similarity  parameter  X  can  be  written 

TL.  ® 

f<%)  r  *  f(T))d71  +  2  r  f3(Tj)dT| 

x  -  ° _ _n* _ 

TL  ® 

f(V  JG  f(7))d71  +  J*  f8(7))dT) 


(1.31) 
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which  is  clearly  a  function  of  the  point  of  zero  shear  stress  T|*.  Glauert 
has  assumed  in  his  analysis  that  the  shear  stress  is  zero  at  the  maximum 
velocity  point.  However,  the  validity  of  this  assumption  has  been  ques¬ 
tioned  by  several  recent  i nvest i gator s  (2,  3,  4)  who  have  found  that  the 
point  of  zero  shear  stress  occurred,  for  the  most  part,  on  the  wall  side 
of  the  maximum  velocity  point  where  the  velocity  gradient  was  largest 
(i.e.,  T|#  <  Tl,,,).  If  it  is  assumed  that  the  point  of  zero  shear  stress 
coincides  with  the  maximum  velocity  point,  then  equation  (1.31),  together 
with  the  assumed  form  for  the  velocity  distribution,  yields 

X  =  1.7977  (1.32) 

The  effect  of  a  different  location  for  the  zero  shear  point 
on  the  similarity  parameter  is  shown  in  Figure  I.  It  can  be  seen  that 
the  value  of  X  is  not  greatly  affected  and  varies  almost  linearly  with 
T)#  provided  the  zero  shear  point  T]^  is  reasonably  close  to  T|m .  The  loca¬ 
tion  of  the  zero  shear  point  wi  II,  however,  have  a  much  stronger  effect 
on  S '  (T|)  and  S(T|)  as  given  by  equations  (1.28)  and  (1.29). 

The  values  of  SC?))  and  S'(T))  corresponding  to  X  =  1.7977  (i.e., 
T|#  =  T|m)  are  tabulated  in  Table  I  along  with  the  similarity  function  f(T|). 
The  effect  of  a  different  location  for  the  zero  shear  point  on  these  func¬ 
tions  can  then  be  obtained  from  the  Taylor  series  expansions 

SU);  X)  =  SUl;  '-7977)  +  jiS.  (T);  1-7977)  (X-  1.7977)  +•••  (1.33a) 

dX 

and 

S’(T);  \)  =  S '  (Tl;  1.7977)  +  dSl  (Tl;  1.7977)  (X  -  1.7977)  +  ••• 

dX  d.33b) 

where  the  derivatives  dnS  and  d’V can  be  obtained  from  equations  (1.28) 

dXn  dXn 


and  (1.29). 


The  characteristic  jet  thickness  6  as  well  as  the  maximum  velo¬ 
city  in  the  jet  profile  Umax  are  functions  of  the  jet  outlet  velocity,  the 
downstream  coordinate,  the  density  and  viscosity  of  the  fluid,  and  the 
height  of  the  jet  outlet.  A  dimensional  analysis  yields 

£  -  fx  (  *,  > 

d  d  v 

or 

6  =  Fj,  (x,  Re0)  (1.34a) 

and 

umax  =  Um  =  Fa  (x,  Re0)  (1.34b) 

Vj 

where  Re0  is  the  jet  outlet  Reynolds  number. 

The  variation  of  these  velocity  and  length  scales  Um  and  5  with 

the  downstream  coordinate  is  partially  governed  by  equation  (1.23).  This 

equation  may  be  integrated  directly  with  the  result 

_  -  I A 

Um  *  (constant)  •  6  (1.35) 

where  X  has  been  assumed  constant.  If  the  region  of  simi larity  in  the 
flow  extended  as  far  upstream  as  the  jet  outlet,  then  the  r-nstant  appear¬ 
ing  in  this  equation  could  be  evaluated  from  conditions  at  the  jet  outlet. 
That  is,  at  x  =0 

Um  =  Ufflfla.  =  I 

Vj 

and 

8=4.=  ! 

d 

so  that  the  constant  would  be  unity.  However,  since  conditions  at  the 
jet  outlet  are  somewhat  dependent  on  the  local  geometry  of  the  outlet 
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as  well  as  the  jet  outlet  Reynolds  number,  an  additional  factor  M  may 
be  introduced  which  takes  into  account  the  effects  of  geometry  and  jet 


Reynolds  number.  Equation  (1.35)  may  then  be  written 


Um  =  ( I  +  M)  6 


,  -l/X 


(1.36) 


From  equation  (1.17),  the  wall  shear  stress  can  then  be  written 

uT*  =  -  o( I  +  M)s  (  6  (X"2)/X  )  (1.37) 

dx 


In  the  case  of  a  two-dimensional  free  turbulent  jet,  it  can  be 
shown  from  momentum  considerations  that  the  maximum  velocity  in  the  jet 
profile  varies  as  x-1'3  while  the  jet  thickness  varies  linearly  with  x, 
the  shear  stress  being  zero  on  the  center  line  by  reasons  of  symmetry. 

An  examination  of  equation  (1.31)  shows  that  as  the  point  of  zero  shear 
stress  T]#  approaches  zero,  X  approaches  a  value  of  two.  Consequently, 
equations  (1.36)  and  (1.37)  both  tend  toward  those  governing  a  free  tur¬ 
bulent  jet  flow. 

Because  the  wall  boundary  layer  represents  only  a  small  portion 
of  the  flow,  the  spreading  and  entrainment  action  of  a  wa I  I  jet  is  quali¬ 
tatively  simi lar  to  that  of  a  free  turbulent  jet.  The  general  forms  of 
the  velocity  distribution  in  the  turbulent  wall  jet  and  turbulent  free 


jet  are  compared  below. 


y 
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While  there  is  no  change  in  momentum  in  the  free  jet  case,  the  wall  jet 
does  suffer  a  momentum  loss  due  to  the  friction  at  the  wall.  The  varia¬ 
tion  of  this  surface  friction  is  gi ven  by  equation  (1.17). 

In  both  cases,  the  total  kinetic  energy  contained  in  the  flow 
up  to  any  point  x  is  given  by 

KE  (x)  -  f  [  f  )dy  ]  dx  (1.38a) 

°  2  3x  2  5y 

Making  use  of  the  fact  that  similarity  also  exists  for  the  free  turbulent 

jet  flow,  this  expression  may  be  integrated  with  the  result 

KE  (x)  *1  p  V,3  d  5  Um3  J  f3  (Tl)dTl  (1.38b) 

2  J  o 

The  amount  of  energy  required  to  insert  the  "wall"  Into  the 
free  jet  flow  is  given  by  J  tw  [  Umax  ]f.EE  dx  where  [  Umax  ]f*EE  is 

JET  JET 

the  maximum  velocity  which  existed  in  the  free  jet  velocity  profile.  It 
is  therefore  possible  to  write 

x 

CKx  6  Um3  ]Epee  “  J*  urS  [Um^EREE  =  Kg  5  Um3  (1.39) 

JIT  J  E  T 

where  the  Kj  (i  ■  I,  2,  •••)  are  constants.  Differentiating  with  respect 
to  x. 


Ki  [  6  Um3  ]p*EE  “  uTa  [Un]ff|EE  =  Kj  ji_  (6  Um3) 
dx  jet  jet  dx 


(1.40) 


JET  JET 

Making  use  of  equations  (1.36)  and  (1.37)  and  the  known  variations  of 

6  and  Um  for  a  free  turbulent  jet,  this  becomes 

„  --3/S  __1/3  _a-2)/X  T(X-3)/X  .  . 

Kax'  +  x  '  d  (6  )  +  1C  d_  (  8  )-0 


SL.  (  6 
dx 


)  +  K*  sL  C  5 

dx 


or 


--2/X  -  „  .-1/3  -"l/*  ,  .-3/3 

J  di[x  +  Ks  8  ]  *  Kj  x 

dx 


(1.41) 
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which  may  also  be  written 


-(X  2)/X  ^  ^  *-i/»  +  5-l/x  ]  =  Kg  j^-a/a  (1.42) 

5  dx 

Previous  experimental  investigations  have  shown  that  5  varies 
almost  linearly  with  the  downstream  coordinate  x  in  the  case  of  a  wall 
jet.  In  order  to  obtain  an  approximate  solution  to  equation  (1.42),  the 
factors  6^  an(j  5  may  be  assumed  to  vary  as  and  x 

respectively.  Then 

_L  Jdi  =  K-  g<*-*)/(2X) 

*  d*  x'l/a  +  Ka  x"l/X 

.  .  ...  ..  _(2-X)/(2X)  .  .  . 

Substituti  ng  z  =x  and  integrating, 

z 

In  6  -  Ke  J  -Z3dz  +  K10 

z  +  K0 

=  K9  [  i3  -  K,z  -  3  K|  +  Kj  In  (z  +  K.)  ]  +  K10 
2  2 


1  7  1/  r  i  —  (2— X)  /X  v  — 
In  6  *  Ks  [  i  x  -  Kt  x 

2 


I  n  terms  of  x, 

+K|  In  ♦  K.)  ] 

+  Kjj  (1.43) 

As  a  further  approximation,  the  second  and  third  terms  in  the  brackets 

may  be  dropped  so  that 

.  jr  _  „_(2-X)/X  „ 

I n  6  —  ox  +  Kx ^ 


or 

6=0  exp  [  Cx^  X^X  ]  (1.44) 

The  constants  C  and  D  appearing  in  this  equation  wi  I  I  depend  on  the  jet 
outlet  Reynolds  number  as  well  as  the  local  conditions  at  the  jet  outlet. 
Referring  back  to  equation  (I.3&)  it  can  be  seen  that  there  are  three 
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factors  (i.e.,  C,  D,  and  M)  describing  the  variation  of  the  velocity  and 
length  scales  which  must  be  determined  experimentally. 

From  equations  (1.36)  and  (1.44),  the  maximum  velocity  decay 
i s  given  by 

Um  =  (I  +  M)  D"l/X  exp  [  -  £  x(2_X)/X  ]  (1.45) 

X 


Likewise,  from  equations  (1.37)  and  (1.44)  the  wall  shear  stress  varia¬ 
tion  can  be  expressed  as 


ux2  =  Tw  -  oC  (I  +  M)2  d<X_2)/X(  X-2  )3  x 


-(2-2X) /X 


exp 


p  vr 


[2^c,(2-X)/x  ] 


(1.46) 


Because  of  the  approximations  involved  in  obtaining  a  solution 
of  equation  (1.41),  the  above  expressions  for  6  and  Um  do  not  tend  toward 
the  corresponding  expressions  for  the  free  turbulent  jet  as  X  approaches 
a  value  of  two.  It  i s  a  I  so  interesting  to  note  that  the  preceding  analy¬ 
sis  does  not  involve  the  idea  of  a  virtual,  origin  or  hypothetical  point 
where  the  jet  is  assumed  to  originate.  This  is  advantageous  because  in 
many  cases  it  has  been  found  that  the  virtual  origin  location  for  the 
maximum  velocity  decay  did  not  coincide  with  that  for  the  jet  thickness 
variation. 


F.  Comparison  With  Glauert's  Theory 

As  mentioned  previously,  Glauert  has  obtained  a  theoretical 
solution  for  the  velocity  distribution  as  well  as  the  variation  of  the 
velocity  and  length  scales  in  a  plane  turbulent  wall  jet.  Glauert's  so¬ 
lution  for  several  values  of  the  parameter  a  are  shown  in  Figure  2  along 
with  the  profile  defined  by  equation  (1.27)  for  purposes  of  comparison. 
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Glauert's  profiles  have  been  constructed  from  the  figures  given  in 
Reference  I  and  no  great  accuracy  is  claimed.  Near  the  wall,  equation 
(1.27)  coincides  almost  exactly  with  the  profile  corresponding  to  a  = 

1.2,  while  far  from  the  wall  it  lies  somewhat  above  Glauert's  solutions 
for  those  values  of  a  which  are  shown. 

The  rate  of  growth  of  the  jet  and  rate  of  decay  of  maximum  velo¬ 
city  are  predicted  by  Glauert's  analysis  as  functions  of  the  Reynolds 
number  based  on  jet  thickness  and  maximum  velocity  in  the  jet  profile. 

If  Umax  and  6  are  assumed  to  vary  as  xa  and  xb  respectively,  then  the 
exponents  a  and  b  are  related  to  the  parameter  a,  which  depends  on  the 
Reynolds  number,  in  two  possible  ways.  That  is, 

a  =  ..-fa _ 

5  +  fa 

b  =  4  +  fa 
5  +  4a 


This  last  set  of  relations  appears  to  be  more  acceptable  ex¬ 
perimentally.  If  the  shear  stress  at  the  point  of  maximum  velocity  were 
zero  as  Glauert  has  assumed,  then  these  two  relations  along  with  equations 
(1.32)  and  (1.35)  imply  that  a  ss  1.25.  Glauert's  profi  le  for  a  =  1.25 
is  shown  in  Figure  3  along  with  equation  (1.27),  and  it  can  be  seen  that 
the  differences  between  the  two  are  not  appreciable. 
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Part  2:  Turbulent  Wall  Jet  Beneath  a  Uniform  Free  Stream 
A.  Similar  Solution 

For  the  more  practical  case  of  a  turbulent  wall  jet  beneath  a 
secondary  uniform  stream,  the  boundary  layer  equations  (l.la)  and  ( I .  lb) 
of  the  preceding  part  are  still  applicable.  The  general  form  of  the  mean 
velocity  distribution  for  this  case  is  as  shown  below,  the  secondary  free 

stream  velocity 
being  denoted  by 

u„. 


Because  the  ratio  of  to  Umax  varies  continuously  in  the  down¬ 
stream  direction  from  a  large  value  near  the  jet  outlet  to  unity  when  the 
jet  disappears,  a  similar  solution  of  the  form  (1.5)  cannot  be  found.  How¬ 
ever,  a  different  type  of  similarity  can  be  found  which  utilizes  the  two 
similarity  functions  found  in  the  analysis  of  the  plane  turbulent  wall  jet. 
That  is,  a  solution  of  the  form 

u  =  Um  f (71)  +  C  gttP  +  h<Tp  +£  j (7))  +  •••  (2.1) 

Um  (4 


may  be  sought  where 


(2.2) 
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and 

T|  =  y/6,  (2.3) 

6  again  being  some  characteristic  thickness  of  the  jet  profile. 

Several  properties  of  the  similarity  functions  g(7|),  h  (T|> ,  j  (T|) ,  *  * 
can  be  immediately  stated,  namely 


g(0)  =  h  (0)  =  J(0)  =  •  •  •  =0 

(2.4a) 

g(“)  =  1 

(2.4b) 

and 

h(®)  =  j  (oo)  =  . . .  =0 

(2.4c) 

At  the  point 

of  maximum  velocity  7]m, 

gCT)m)  -  h(7]m)  =  J(V  *  •••  =0 

(2.5a) 

and 

g'(Tlm)  =  h*  (Tin,)  -  J' <1)m  -  •••  =0 

(2.5b) 

Integration  of  the  continuity  equation,  making  use  of  equation 
(2.1),  resu Its  in 

7] 

v  =  Um  di  I]  f<7])  -  jL  (5Um)  J  f (z)dz  +  C  dj.  [  Tj  gCT)) 
dx  dx  0  dx 

-  f  9(z)dz  ]  +  Ca  [  1_  d£  7]  h  (71)  -  d_  (6/Um)  f  h  (z)dz  ] 

0  Um  dx  dx  “ 

+  C3  [I  di.T)  j  (T))  “  d_  (5/Um3)  J11  j(z)dz]  +  •••  (2.6) 

Uma  dx  dx  ° 

B.  Shear  Simi laritv 

With  the  aid  of  equations  (2.1)  and  (2.6),  the  momentum  equation 


may  be  written 
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Um  fa(11)  -  Um  Jl_  (Sum)  f  ’  (T))  J  f(z)dz  +  C  C  Ba  f(1p  g(T|) 
dX  {  dx  o  dx 

T]  _  m 

-  Um  d£  f '  (Tl)  f  g(z)dz  -  I  i.  (6Um) g '  (T|)  f  f(z)dz  ] 

5  dx  o  6  dx  410 


-  cs  [  l  9 '  C-Tl)  f1  g(z)  dz  +  _J_  £_  (&W  h'(T|)  J1  f(z)dz 

6  dx  °  Um5  dx  o 


+  Urn  i_  (5/Um)  f  '(Tp  P  h(z)dz  ]  -  C3  /  _L_  j^m  [  f  (7p  j  (T)) 
5  dx  do  l  U_a  dx 


Ua  dx 


+  g(T))h(T))]  +11  C8Um>  j '  CTl)  J  f(z)dz  +  I  j&h'CTl) 
SUm3  dx  o  5llm  dx 

Tl  T! 

•  J  g(z)dz  +  ±  (§/Um)  g'  (Tl)  J  h(z)dz  +  Um  sL.  (6/Um3) 

o  5  dx  o  5  dx 


f  ’  (Tl)  f  J(z)dz}  + 


2H*1 

ay 


(2.7) 


Assuming  that  the  shear  stress  in  the  free  stream  is  zero,  this  equation 
may  be  integrated  with  respect  to  y  from  zero  to  ®  i n  order  to  determine 
the  shear  stress  at  the  wall.  That  is, 

uTa  »-al  (5Uma)  -  C  i_  C5Um)  [  2  f  f  (Tp  gCTpdTl  -  j"  f  (T))  dTl  ] 

_»_T  jt.  w  ft  w  - 


dx 


dx 


-  ca  dfi.  [  J*  gs  (if)  dTl  -  r  g(Tl)  dTl  +  2  J*  f  (Tl)  h  <Tl>  dTl  ] 

dx  o  o  o 

-  C3  JL  (5/Um)  C  2  f  f  (Tl)  j  (Tl)  dT)  +  2  f  g  (Tl)  h  CTJ)  dTl 

dx  o  o 

-  f  h (T))  dT)  ]  -  ••• 


(2.8) 


Def I ni ng 


P  ■  Jo  f  (Tl)  C  2  g(Tj)  -  I  ]df) 

y  *  f  g (Tl)  [  g(Tl)  -  l  ]«m  +  2  f  f  (71)  h  Op  dTl 
o  o 

a  ■  £°h(T|)  [  2g(Tl)  -  I  ]dT|  +  2  J*"  f  (Tl)  J(T|>dT| 


(2.9) 

(2. 10) 

(2.11) 
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then  equation  (2.8)  may  be  written 

V  =  '“4  (5lJmS)  -  C  3  (SUm)  -  CaY  4 
dx  dx  dx 

-  C3o  (S/Um)  -  ••  (2. 12) 

dx 

Likewise,  the  skin  friction  coefficient  defined  by  equation  (1.18)  may 
be  written 

cf  =  -  2sk  jL  (5Um3)  -  24  SL  <5um)  -  2£!v  jfl, 

Uma  dx  Um  dx  Ums  dx 

-  l  i_  (5/Um)  -  •••  (2.  13) 

Um3  dx 

Assuming  that  similarity  exists  for  the  shear  stress  in  a  wall 

jet  beneath  an  external  free  stream,  then  analogous  to  equation  (2.1) 

u„2  *  ur3  CSOp  +£_Tnp  +£i  W (T|)  +j£Z(7p  +  •••  (2.14) 

Um  Uma  U/ 

This  equation  may  then  be  combined  with  the  momentum  equation  (2.7).  |n 
order  that  the  resulting  equation  be  valid  for  all  values  of  £,  the  coef¬ 
ficient  of  each  £n  must  vanish.  For  n  =  0 

Um  dUm  [  f2  (T))  -  f  *  (71)  f(z)dz  +2aS'(Tl)  ]  = 

dx  o 

V  d|  [  f  ’  (Tp  J71  f(z)dz  -  aS'(T|)  ]  (2.15) 

?  dx  o 

which  is  the  same  as  equation  (1.22). 


Tl 

6  dUm  f  (Tp  g  (Tj)  -  Um  4  f'Clp  J  g(z)dz  -  (6Um)g*  (7]) 

dx  dx  °  dx 

•  f(z)dz  =  -  8  (6Um)S '  (T|)  -  a_iL  (6U*>T'  (71)  (2.16) 
°  dx  Um  dx 


For  n  *  I 
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For  n  =  2, 


AL  g'(Tl)  J  g(z)dz  +  J_4_  (3um)  h '  (T|)  J*11  f(z)dz  +  Um  (5/Um) 


dx 


Um  dx 


dx 


f'(Tl)  f  1  h(z)dz  =  y  AL  S' (TP  +  &_  A_  (6Um)  T*  (Tp 
o  dx  ■■  ■" 


Um  dx 


+  _flL  A-  (SUm2)  W'  <T1) 

1 1  2  j- 


(2. 17) 


Um  dx 


For  n  =  3, 


J\ 


_6  dUm  [  f  (T|)  j  (1))  +  g(7])h(T])  ]  + _ I  (5Um)  j  '  (T))  f  f(z)dz 

Um*  d*  "  2  ~ 


Um"  dx 


J) 


+  J_  45.  h '  (7))  J”  g(z)dz  +  £_  (6/Um)g'nj)  f  h(z)dz  +  Um 

Um  dx  o  dx  o  dx 

•  (6/Um2)  f'(Tl)  f1  j(z)dz  =  a  ±_  (6/Um)  S’  (7))  +  T '  (T)) 

o  dx  Um  dx 

+  _fi Ti  <6Um)  W*  (71)  +  _fld  (5Um2)  Z 1  (7))  (2.18) 

Um  dx  V  dx 

The  boundary  conditions  governing  the  shear  similarity  functions 
T  (T|) »  W(T|),  Z(7p,  •••  may  be  written 

T (0)  =  W(0)  =  Z(0)  =  • • •  =  0  (2.19a) 

T(oo)  =  w(co)  =  z(«)  =  •  •  •  =  0  (2.19b) 

T'(0)  =  W'(0)  =  Z '  (0)  =  •••  =  0  (2.19c) 

In  order  that  there  be  a  similarity  solution  independent  of  the 
downstream  coordinate,  each  of  the  equations  (2.15)  through  (2.18)  must 
be  independent  of  x.  As  before,  equation  (2.15)  implies  the  condition 

dUm  +  J.  =  0  (2.20) 

Um  dK  5  dx 

with  solution  given  by  equation  (1.35).  Equations  (2.15)  through  (2.18) 
may  then  be  written 

fs <T|)  +  (X-l)  f » (71)  J11  f  (z ) dz  +  (2-\)  a  S '  (71)  =  0 


(2.21) 
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fOp  g(T|)  +X  f  *  <T1)  g(z)dz  +  (X- 1 )  g'(Tl)  f  Uz)dz 

o  o 

+  ( I  -X ) e  S'  (11)  +  (2-X)a  T'(Tp  =  0  (2.22) 

Xg'  (T|)  f11  g(z)  dz  +  (X-l)  h  •  (-fl)  rT|  f(z)dz  +  (X  +  D  f'CTp 
"o  *o 

Jl) 

•  J  h(z)dz  -  Xy  S '  (71)  +  (l-X)g  T'(T)) 

o 

+  (2-X)a  W'(Tp  =  0  (2.23) 

T|  T| 

f(Tpj(71)  +  g(T|)h(Tl)  +  (l-X)j'(lp  J  1  f(z)dz  -  Xh '  (Tp  J  g(z)dz 

O  0 

-  (X  +  D  g'  (71)  f  h(z)dz  -  (2+X)  f '  (T|)  f1  j (z)dz  +  (X+I)ct 
o  o 

•  S'(Tp  -  XyT'(Tl)  +  (I-X)p  W'(T1)  +  (2-X)aZ'(T))  =0  (2.24) 

These  are  the  similarity  equations  which  govern  the  flow  in  a  wall  jet 
beneath  a  uniform  free  stream. 

If  the  point  of  zero  shear  stress  is  again  denoted  by  T]#,  then 
in  addition  to  equations  (2. 19) 

T (71)  =  W(T1  )  =  Z (T|  )  =  ■••  =0  (2.25) 

w  IT  IP 

It  might  a ! so  fee  pointed  out  that  equation  (1.25)  relating  the  simi larity 
parameter  X  to  T|*  is  sti  I  I  valid  in  the  case  of  a  wa  I  I  jet  beneath  a  secon¬ 
dary  free  stream. 

C.  Simi laritv  Functions 

Elimination  of  S '  <T|)  from  equations  (2.21)  and  (2.22)  gives 

f(T])  g(T|)  +  X  f '  (7])  f11  g(z)dz  +  (X-l)  g'(ip  J11  f(z)dz  +  (l-X)B 

o  o  (X-2)a 

T] 

•  Cfs(Tp  +  (X-l)  f  ’  (*n)  J  f(z)dz]  +  (2-X)aT'(Tp  =  0 


(2.26) 
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As  before,  in  order  to  obtain  a  solution  of  this  equation,  some  further 
assumptions  must  be  made  concerning  the  velocity  distribution  g(T|)  or  the 
shear  distribution  T (Tp .  One  possible  form  for  the  function  g(7p  which 
satisfies  the  boundary  conditions  (2.4)  and  (2.5)  and  which  is  suggested 
by  the  function  f(Tp  found  in  the  preceding  part  is 

g(T|)  =  {erf(BT))  [  I  -  e"1^  ]4  [  I  -  f(7|)]}2  (2.27) 

where  B  has  the  value  .67753  as  determined  in  the  preceding  part. 

From  equation  (2.26) 

T| 

T'  (Tp  =  (l-X)B  [  f2 (7))  +  (X-D  f  ’  (71)  J*  f  (z)dz  ]  + _ I _ 

[ (X-2)a]a  o  (X-2)a 

7]  m 

•  [f  (7|)g(7])  +Xf'(T))  J*  g(z)dz  +  (X- 1  )g*  (Tp  J  f(z)dz]  (2.28) 

o  o 

which  may  be  integrated  giving 

T (T|)  ■  (I -X) B  [ (2-X)  f71  f3(z(dz  +  (X-I)f(7p  f ^  f(z)dz 
[0-2)0?  o  o 

+  I  [2 <  I -X )  f  f(z)g(z)dz  +  Xf(T|)  g(z)dz 

(X-2)a  o  o 

+  (X- I )  g(Tp  ^  f (z) dz]  (2.29) 

o 

If  S '  (TJ)  and  T '  (71)  are  eliminated  from  equation  (2.23),  there 

resu Its 

Xg' (71)  J71  g(z)dz  +  (X-Dh'(Tl)  J71  f(z)dz  +  (X+i)f'(7))  J  h(z)dz 
o  o  0 

+  f  r  (X-l)Bl2  -  Xv  )  [  f 3  ( T|)  +  (X- 1  j  f  •  (7|>  r71  f(z)dz] 

U(x-2)aJ  (X-2)aJ  0 

+  (l-X)B  [f (Tj)g(Tl)  +  Xf'(71)  g(z)dz  +  (X-l)g'(Tl)  f71  f(z)dz] 

(X-2)a  o 


+  (2-X)a  W'  (7j)  -  0 


(2.30) 
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A  reasonable  form  for  the  function  h(7|)  satisfying  the  required  boundary 
conditions  (2.4)  and  (2.5)  is 


h ( T|)  =  50  f (T))  g (T|)  (2.31) 


From  equation  (2.30) 

11  v 

WO))  =  __l_  [Xg'(T)>  I*  g(z)dz  +  (X-l)  h '  (71)  J  f(z)dz  +  (X+Of'CTp 
(X-2)a  o  o 

•  /  h(z)dz1  +  (l-X)B  [f(T|)g(T])  +  Xf'(Tl)  f1  g(z)dz  +  (X-l) 

o  C(X-2)a?  o 

J| 

•  g'(Tl)  S  f(z)dz]  + _ L_^  f  r(X-l)8ls  -  X  V  )  rf8(Ti) 

°  [ (X-2)a]a  L  (X-2)a  / 

T\ 

+  (X-l)  f '  (71)  J  f(z)dz]  (2.32) 


which  may  be  integrated  with  the  result 

T|  v  Tj 

W (T|)  = _ I _  [Xgcrp  f  g(z)dz  -  X  J  ga(z)dz  +  (X-l)h(Tl)  J1  f(z)dz 

(X-2)a  o  o  o 


+  (X  +  l)  f (71)  f1  h(z)dz  -  2X  f71  f(z)h(z)dz]  +  (l-X)B 

Jo  [ (X-2)a]3 

•  [2(  I  -X)  J71  f  (z)g(z)dz  +Xf(T])  J11  g(z)dz  +  (X- 1 )  g  (T))  f  f(z)dz] 
o  n  ° 


T) 


+ _ ! _ ^  [  r(X-l)Bl8  -  X  v  \  r (2-X)  r  faU)dz 

[(X-2)aJs  l  (X-2)a  /  o 


+  (X-l)  f (Tp  J71  f(z)dz] 


C2.33) 


o 

Following  a  procedure  similar  to  the  above,  S'(T)),  T '  (Tj) ,  and 
W'(T])  may  be  eliminated  from  equation  (2.24)  by  means  of  equations  (1.28), 
(2.28),  and  (2.32).  The  resulting  equation,  however,  still  contains  the 
two  unknowns  j  (T|)  and  Z '  (TJ)  and  because  of  its  length  will  not  be  stated 
here.  One  possible  form  for  the  function  j(T|)  which  satisfies  the  pre¬ 
scribed  boundary  conditions  is 
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J (Ip  =  160  e'BT,S  f <T|)  h (T|)  (2.34) 

With  this  assumption,  the  function  Z (T))  may  then  be  obtained  by  a  straight¬ 
forward  integration. 

An  examination  of  equation  (2.1)  shows  that  a  series  representa¬ 
tion  of  this  form  is  valid  only  up  to  the  point  where  the  Jet  disappears 
(i.e.,  only  when  (£/Um)  =  <U00/Umax)  <  I).  As  long  as  £/Um  remains  less 
than  one,  this  series  wi I  I  converge  provided  each  of  the  simi larity  func¬ 
tions  is  bounded.  Beyond  the  point  where  the  jet  disappears,  the  flow 
becomes  equivalent  to  a  normal  turbulent  boundary  layer  in  a  zero  pressure 
gradient. 

Since  the  normal  region  of  interest  for  the  wall  jet  flow  beneath 
an  external  free  stream  extends  only  up  to  a  point  where  the  ratio  C/Um 
is  still  somewhat  less  than  unity,  the  velocity  distribution  can  be  ade¬ 
quately  represented  by  the  first  four  terms  of  equation  (2.1).  That  is, 

!i_  =  f  (T|)  +£_  g(Tp  +£?.  h (T|)  +£L  J(Tp  (2.35) 

Urn  Um  Um2  Um3 

Likewise,  the  shear  stress  distribution  can  be  adequately  represented  by 

u®  =  u_a[S(Tl)  +  £_  T(7))  +  II  W(T])  +  Z (T|) ]  (2.36) 

Um  Um3  Um3 

The  next  term  in  each  of  these  expressions  becomes  important  only  in  the 
case  where  £  is  relatively  large,  or  very  far  downstream  of  the  jet  out¬ 
let  where  Um  becomes  small. 

The  simi  larity  functions  g (T|) »  h(T|),  and  j(T|)  were  programmed 
for  the  computer  along  with  the  function  f(T))  of  the  preceding  part. 

Values  of  these  functions  are  listed  in  Tables  El,  HI,  and  DT.  The 
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constants  (3,  V/  and  &  defined  by  equations  (2.9),  (2.10),  and  (2.11)  were 
found  to  be 

00 

0  =  J*  f  (Tj)  [2g(Tl)  - 1  ]dT|  =  '  i-oliz  (2.37) 

o 

Y  =  [g(T|>  -l]dTl  +2  f  f(T))h(71)d71  =  -.0336  1  (2.38) 

”  o  0  o 

and 

00  00 

a  *  J*  h(Tl)  [2g(T))  -l]dT]  +2  J*  f (T|)  j (T|) dT)  =  +.63zro  (2.39) 

o  o 

The  shear  simi  larity  functions  T (T|) ,  WOT]),  and  Z(T))  together 
with  their  derivatives  are  also  listed  in  the  above  tables.  These  values 
were  computed  under  the  assumption  that  the  point  of  zero  shear  stress 
coincides  with  the  maximum  velocity  location.  The  effect  of  a  different 
location  for  the  zero  shear  point  on  these  functions  (i.e.,  \  £  1.7977) 
can  be  obtained  by  simple  Taylor  series  expansions  as  outlined  in  the  pre¬ 
ceding  part. 

D.  Variation  of  Velocity  and  Length  Scales 

The  presence  of  an  external  free  stream  above  the  wall  jet  alters 
the  functional  dependence  of  6  and  Um  as  given  by  equations  (1.34a)  and 
(1.34b).  To  account  for  the  secondary  stream,  these  equations  may  be  re¬ 
written  in  the  form 

6  =  Fx (x,  £,  Reffl )  (2.40a) 

and 

Um  =  F8(x,  Q,  Re0 )  (2.40b) 

As  was  pointed  out  earlier,  equation  (1.35)  still  relates  the  velocity 
and  length  scales  even  with  the  addition  of  the  secondary  free  stream. 
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Due  to  the  greater  momentum  in  the  jet  profl le,  mass  is  entrained 
into  the  Jet  from  the  free  stream  but  not  to  as  great  an  extent  as  when 
the  free  stream  is  not  present.  Thus  a  wall  jet  necessarily  mixes  more 
slowly  with  a  moving  stream  than  with  still  air  and  only  for  very  small 
values  of  £  will  the  rate  of  growth  of  the  jet  approach  that  of  the  plane 
turbulent  wall  jet.  Analogous  to  the  stl  II  air  case,  the  most  obvious 
choice  for  the  characteristic  jet  thickness,  6,  is  that  value  of  y  in  the 
outer  portion  of  the  jet  where  the  velocity  lies  halfway  between  the  maxi¬ 
mum  and  free  stream  velocities.  This  point  is  also  measured  very  easi  ly 
exper imenta lly. 

An  examination  of  Tables  I  through  II  reveals  that  the  simi larity 
functions  as  tabulated  do  not  take  into  account  the  fact  that  the  jet  thick¬ 
ness  is  also  dependent  on  the  velocity  ratio  £.  That,  y  =  6  does  not  nec- 
essari  ly  correspond  to  T)  =  I  except  for  the  case  £  =0  which  is  that  of 
the  plane  turbulent  wall  Jet.  This  fact,  however,  can  be  used  to  deter¬ 
mine  the  effect  of  the  external  free  stream  on  the  velocity  and  length 
scales.  Denoting  the  value  of  7)  in  equations  (1.26),  (2.27),  (2.31),  and 
(2.34)  where  the  velocity  lies  halfway  between  the  maximum  and  free  stream 
velocities  by  \/z,  then  from  equation  (2.35) 

id  +C  )  =  f(Tll/a)  +  5  g(%/8)  +  5l  h(1i/3)  +51  jHll/a)  (2.41) 

2  Um  Um  Um2  Um 

This  relation  between  \ and  £/Um  is  shown  in  Figure  4  and  can  be  written 

symbolically  in  the  form 

\/2  =  F(£/Um)  (2.42) 

In  constructing  a  given  velocity  profile  by  means  of  the  Tables,  each  value 
of  T|  must,,  therefore,  be  scaled  by  the  factor  [F(£/Um)l_1  appropriate  to 
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the  given  profile.  An  example  velocity  profile  construction  is  given  in 
the  Appendix. 

The  introduction  of  the  scale  factor,  F,  necessarily  modifies 

equation  (2.20)  which  must  now  be  written 

L_  fj^m  +  j _  d_  (5F)  =  0 

Um  dx  5F  dx 

with  solution 

Um  =  const.  (6F)"l/X  (2.43) 

Analogous  to  equation  (1.36)  for  the  plane  turbulent  wall  jet,  this  may 
be  rewritten  in  the  form 

Um  =  (UM)  (6F)‘l/X  (2.44) 

where  the  factor  M  again  takes  into  account  local  conditions  at  the  jet 

outlet.  Likewise,  the  wall  shear  stress  must  be  expressed  as 

uTa  =  -  a  i.  (6FUm3)  -  CP  JL  (6FUm)  -  Cs  y  d_  (6F> 
dx  dx  dx 

-  C3  a  £_  <££)  -  •••  (2.45) 

dx  Um 

It  may  be  pointed  out  that  these  results  are  consistent  with  the  results 
for  the  plane  turbulent  wall  jet,  since  the  scale  factor  F  is  unity  for 

C  =  0. 

The  fact  that  the  values  of  T]  i n  the  Tables  are  correct  only  for 
C  =  0  and  must  be  scaled  by  the  factor  F  1  for  £  £  0  implies  that  the  ef¬ 
fect  of  the  external  free  stream  on  the  growth  of  the  jet  is  given  by  the 
function  F(C/Um).  From  equation  (1.44),  the  jet  thickness  variation  is 
therefore  given  by 

fi  =  D  F<C/Um)  exp  [Cx(2"X,)/X]  (2.46) 


I 
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while  the  maximum  velocity  decay  is  given  by 

Um  -  (I  +  M)  [D  Fs(C/Um)]‘l/X  exp[  -  £  x(2_X)/X  ]  (2.47) 

X 

Again,  the  constants  C  and  D  appearing  in  these  equations  wi II  depend  on 
conditions  in  the  neighborhood  of  the  jet  outlet. 

If  the  above  equations  (2.46)  and  (2.47)  are  rewritten  in  the 

form 

i  =  D  exp[C  x(2_X)/X  ]  (2.48) 

F 

and 

Um  Fa/X  =  (I  +  M)  D-l/X  exp[  -  £  x(2'X)/X  ]  (2.49) 

X 

then  the  quantities  on  the  left  are  dependent  on  the  velocity  ratio  C  only 
to  the  extent  that  the  factors  C,  D,  and  M  are  dependent  on  £.  Except  for 
the  case  Q  =  0,  a  solution  of  these  equations  for  given  values  of  C,  D,  and 
M  must  be  obtained  by  trial  and  error  in  conjunction  with  Figure  4  which 
describes  the  variation  of  F  with  £/Um. 
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131  EXPERIMENTAL  INVESTIGATION 

The  experimental  phase  of  this  investigation  was  carried  out 
in  the  Princeton  University  2"  x  36"  two-dimensional  smoke  tunnel  which 
has  a  maximum  velocity  capability  of  approximately  90  feet  per  second. 

A  portion  of  the  test  section  floor  of  this  tunnel  was  removed  in  order 
to  accommodate  the  jet  and  rear  plate.  The  entire  model  was  machined  from 
aluminum  and  a  sectional  view  of  its  installation  is  shown  in  Figure  5. 

The  rear  plate,  or  wall  as  it  is  commonly  called,  consists  of  two  pieces 
and  has  an  overall  length  of  34.5  inches.  The  height  of  the  jet  outlet 
could  be  varied  between  .03  inches  and  .15  inches  by  merely  raising  or 
lowering  the  rear  plate.  Sealing  between  the  entire  model  and  tunnel  walls 
was  accomplished  by  means  of  "0"  ring  segments  and  rubber  fi I lers. 

Located  upstream  of  the  jet  outlet  is  a  suction  slot  which  pro¬ 
vides  some  control  over  the  boundary  layer  which  exists  above  the  jet  out¬ 
let  when  the  external  free  stream  is  used.  The  present  investigation, 
however,  does  not  include  any  measurements  to  determine  the  effect  of  this 
initial  boundary  layer  on  the  flow  field  downstream  of  the  jet  outlet. 

In  each  case  in  which  the  external  stream  was  used,  a  normal  flat  plate 
boundary  layer  was  allowed  to  form  on  the  jet  lip. 

The  air  issuing  from  the  jet  was  supplied  by  two  centrifugal 
compressors  connected  in  series.  To  avoid  any  angularity  in  the  flow  with 
respect  to  the  plane  of  symmetry  of  the  tunnel,  a  small  honeycomb  element 
was  installed  in  the  plenum  chamber  just  upstream  of  the  jet  outlet.  With 
this  arrangement,  the  velocity  distribution  at  the  jet  outlet  was  found 
to  be  two-dimensional  to  within  +  and  the  maximum  Reynolds  number 
based  on  jet  velocity  and  jet  slot  height  which  could  be  attained  was  in 
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the  neighborhood  of  21,000.  Both  the  plenum  chamber  and  suction  chamber 
were  effectively  sealed  by  means  of  "0"  ring  segments.  The  height  of  the 
jet  outlet  could  be  held  to  within  +  .001  inches  of  any  desired  value  by 
means  of  a  set  of  cylindrical  plug  guages. 

The  initial  phase  of  the  experimental  program  consisted  mainly 
of  a  verification  of  the  fact  that  the  flow  remained  sufficiently  two- 
dimensional  for  a  reasonable  distance  downstream  of  the  jet  outlet.  Trans¬ 
verse  velocity  distributions  in  the  Jet  were  measured  at  various  downstream 
locations  both  with  and  without  the  external  free  stream.  This  was  accom¬ 
plished  with  the  aid  of  a  rake  consisting  of  ten  evenly  spaced  total  head 
tubes  .032  inches  in  diameter,  the  pressures  being  measured  by  means  of 
a  bank  of  alcohol  manometers,  in  all  cases,  it  was  found  that  the  flow 
remained  two-dimensional  to  within  +  5%  in  a  central  lane  approximately 
one  inch  wide  extending  as  far  as  20  inches  downstream.  Consequently, 
a  I  I  subsequent  measurements  were  carried  out  in  this  narrow  region. 

In  order  to  determine  the  static  pressure  variation  in  the  longi¬ 
tudinal  direction,  the  rear  plate  of  the  model  was  equipped  with  a  series 
of  static  pressure  taps.  The  location  at  which  the  longitudinal  pressure 
gradient  essentially  became  zero  could  then  easi  ly  be  determined,  and  all 
measurements  were  performed  well  downstream  of  this  point. 

The  velocity  profiles  normal  to  the  rear  plate  were  measured 
by  means  of  a  total  head  tube  having  an  internal  diameter  of  .020  inches 
and  an  external  diameter  of  .032  inches.  This  probe  was  mounted  on  a  car¬ 
riage  as  shown  in  Figure  6  as  was  the  total  head  rake  used  in  the  two- 
dimensionality  study.  The  carriage  was  guided  during  its  fore  and  aft 
movement  by  a  small  track  which  is  an  integral  part  of  the  rear  plate. 
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The  long  threaded  rod  extending  downstream  from  the  rear  of  the  carriage 
was  used  to  control  the  longitudinal  position  of  the  carriage  and  total 
head  tube.  The  vertical  movement  of  the  probe  was  powered  by  a  small  di¬ 
rect  current  electric  motor  mounted  on  the  carriage,  the  detai Is  of  the 
drive  mechanism  also  being  shown  in  Figure  6.  The  vertical  position  error 
of  this  system  was  found  to  be  less  than  .002  inches  over  the  entire  range 
of  movement  of  the  probe.  When  the  external  free  stream  was  employed,  the 
pitot  static  tube  located  at  the  top  of  the  carriage  was  used  to  measure 
its  velocity. 

Initially,  the  pressures  sensed  by  the  total  head  probe  were 
fed  to  an  inclined  manometer  whose  angle  of  inclination  could  be  varied 
between  7  and  45  degrees.  Data  for  several  velocity  profiles  were  re¬ 
corded  using  this  manometer,  but  it  was  found  that  because  the  length  of 
time  involved  in  making  such  measurements  was  so  long,  any  appreciable 
change  in  ambient  conditions  could  substantially  alter  the  results.  Con¬ 
sequently,  the  data  for  the  remaining  profiles  were  recorded  with  the  aid 
of  a  differential  pressure  transducer  having  a  range  of  +  .6  psi.  Because 
of  the  nearly  instantaneous  response  of  such  a  device,  the  measurement 
time  was  drastically  reduced  without  any  appreciable  loss  in  accuracy. 
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II  DISCUSSION 

A.  Comparison  of  Experimental  Results  With  Theory 

I .  Plane  Turbulent  Wa I  I  Jet 

Ve loci ty  prof i les  in  the  p lane  turbu lent  wa I  I  jet  ( i .e. ,  Q  =  0) 
were  measured  for  two  different  values  of  the  jet  outlet  Reynolds  number. 
These  profiles  are  shown  in  Figures  7  and  8  along  with  the  similarity  func¬ 
tion  f(T|)  defined  by  equation  (1.27)  for  comparison  purposes.  These  re¬ 
sults  apparently  verify  the  assumption  that  similarity  with  respect  to 
the  velocity  profiles  exists  at  least  as  far  as  417  slot  widths  downstream 
of  the  out  let. 

The  data  presented  in  Figure  7  were  recorded  with  the  aid  of  the 
inclined  manometer  discussed  previously  while  those  in  Figure  8  were  re¬ 
corded  by  means  of  the  pressure  transducer.  Extreme  care  was  taken  in 
making  the  measurements  in  the  outer  portion  of  the  jet  where  the  velocity 
tends  to  zero.  There  is  an  appreciable  source  of  error  in  this  region 
due  to  the  fact  that  the  turbulence  level  is  very  high  and,  because  of  the 
mass  entrainment,  the  transverse  inflow  velocity  is  relatively  large  re¬ 
sulting  in  a  yawed  flow  relative  to  the  total  head  tube.  Both  of  these 
effects  combine  to  cause  the  total  head  tube  to  read  low  in  the  outer  poi 
tion  of  the  jet.  There  is  also  a  source  of  measurement  error  in  the  region 
within  a  few  total  head  tube  diameters  from  the  wall.  Here,  errors  may 
arise  from  viscous  effects,  the  fact  that  there  is  a  velocity  gradient 
across  the  face  of  the  tube,  and  the  displacement  effect  due  to  the  pre¬ 
sence  of  the  wall.  Since  these  latter  two  effects  tend  to  cancel  each 
other,  the  overall  error  in  this  region  is  not  appreciable  and  no  correc¬ 
tions  to  the  data  were  employed. 
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As  a  further  verification  of  its  usefulness  in  predicting  the 
turbulent  wall  jet  velocity  profiles,  the  data  of  several  other  investi¬ 
gations  have  been  compared  with  the  similarity  function  defined  by  equa¬ 
tion  (1.27).  The  data  of  George  (7)  and  Myers,  Schauer,  and  Eustis  (8)  are 

presented  in  Figures  9  and  10  respectively  and,  because  of  the  wide  range  of 

downstream  positions  over  which  the  data  were  taken,  they  offer  further  sup¬ 
port  for  the  similar  velocity  profiles  assumption.  The  data  of  Patel  (II), 
Sigalla  (12),  Schwarz  and  Cosart  (2),  and  Forthmann  (10)  covering  a  somewhat 

more  limited  range  of  downstream  positions  are  presented  in  Figures  II,  12, 

13,  and  14  respectively.  Except  for  the  data  of  George,  all  data  were 
taken  from  the  references  cited  and  are  thereby  limited  in  accuracy.  It 
is  easily  seen  that  the  over-all  agreement  of  the  assumed  form  for  the  velo¬ 
city  distribution  with  the  experimental  results  of  this  investigation  as 
well  as  with  those  of  other  previous  investigations  is  remarkably  good. 

Only  in  the  extreme  outer  portion  of  the  profile  is  there  any  appreciable 
discrepancy,  the  data  points  lying  consistently  below  the  curve  as  expected. 

As  a  check  on  equation  (1.44),  the  approximate  solution  found 
for  the  variation  of  jet  thickness  with  downstream  coordinate,  the  data 
have  been  plotted  in  the  form  In  6  versus  x^  ^  \  These  data  are  shown 
in  Figure  15  where  X  =  1.7977  has  been  assumed.  Likewise,  the  data  of 
George,  Myers,  Schauer,  and  Eustis,  Patel,  and  Forthmann  for  x  greater 
than  20  are  presented  in  Figures  16,  17,  and  18  in  the  same  manner.  The 
data  of  Schwarz  and  Cosart  were  not  available  while  the  data  of  Sigalla 
were  not  easily  adaptable  to  the  above  manner  of  presentation.  The  con¬ 
stants  C  and  0  appearing  in  equation  (1.44)  were  then  determined  to  give 
a  reasonable  fit  to  the  data  in  each  case.  The  value  of  (l+M)  =  Um5 * , 
as  defined  by  equation  (1.36),  was  determined  by  simply  averaging  the  values 
of  in  each  case. 
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The  variation  of  $  and  Um  with  the  downstream  coordinate  x  for 
each  of  the  Investigations  mentioned  above  Is  presented  in  Figures  19 
through  26.  It  is  easi ly  seen  from  these  figures  that  6  varies  almost 
linearly  with  x,  which  means  that  the  influence  of  the  wall  is  fairly 
small  since  the  jet  thickness  varies  linearly  with  x  for  the  free  turbu¬ 
lent  jet.  It  Is  also  evident  that  the  assumption  X  =  1.7977  (i.e.,  that 
the  point  of  zero  shear  stress  coincides  with  the  maximum  velocity  loca¬ 
tion)  is  reasonable  as  far  as  the  prediction  of  the  growth  of  the  jet  and 
decay  of  maximum  velocity  is  concerned. 

In  Figure  27,  the  transverse  shear  stress  calculations  of 
FSrthmann  are  compared  with  the  shear  si  mi  larity  function  S (T))  dPfined 
by  equation  (1.29)  for  X  =  1.7977.  Only  the  data  corresponding  to  values 
of  x  greater  than  20  are  shown  since  the  similarity  assumption  is  not  valid 
for  smaller  values  of  x.  The  agreement  between  the  two  is  reasonable  and 
adequately  verifies  the  assumption  of  shear  stress  similarity. 

2.  Turbulent  Wall  Jet  Beneath  a  Uniform  Free  Stream 

Velocity  profi les  in  the  turbulent  wall  jet  beneath  a  secondary 
uniform  stream  were  measured  for  two  different  values  of  the  jet  outlet 
Reynolds  number  which  also  corresponds  to  two  different  values  of  the  velo¬ 
city  ratio  Q.  These  profi  les  are  shown  in  Figures  28  through  32,  along 
with  the  similarity  expansion  defined  by  equation  (2.35)  for  comparison 
purposes.  All  of  these  data  were  recorded  with  the  aid  of  the  pressure 
transducer  discussed  previously.  These  results  verify  that  an  expansion 
in  terms  of  similarity  functions  as  defined  by  equation  (2.35)  Is  certainly 
valid  for  at  least  400  slot  widths  downstream  of  the  outlet  provided  the 
ratio  C/Um  remains  approx i mate  I y  less  than  .4. 
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The  addition  of  a  secondary  stream  considerably  reduces  the  ex¬ 
perimental  errors  arising  in  the  outer  portion  of  the  jet  as  a  result  of 
pitot  tube  incidence.  Consequently,  the  outer  portion  of  the  velocity 
profile  for  a  turbulent  wall  jet  beneath  a  secondary  stream  exhibits  a 
longer  'tail'  than  a  profile  for  a  wall  jet  issuing  into  a  region  of  sta¬ 
tionary  fluid.  However,  the  sources  of  measurement  error  in  the  region 
near  the  wall  still  remain,  but  again  the  over-all  error  in  this  region 
is  not  appreciable  and  can  be  neglected. 

Again,  as  a  further  verification  of  its  usefulness  in  predicting 
the  velocity  profiles  in  the  turbulent  wall  jet  beneath  a  secondary  uni¬ 
form  stream,  the  data  of  several  previous  investigations  have  been  compared 
with  the  velocity  profiles  constructed  on  the  basis  of  equation  (2.35). 

The  measurements  of  Eskinazi(l3)  at  three  different  values  of  jet  outlet 
Reynolds  number  and  velocity  ratio  Q  are  shown  in  Figures  33  through  42. 
These  measurements  were  made  in  a  range  extending  as  far  as  465  slot  widths 
downstream  of  the  jet  outlet.  The  data  of  George  (5),  extending  as  far 
as  400  slot  widths  downstream  of  the  jet  outlet,  are  shown  in  Figures  43, 
44,  and  45.  The  data  of  Patel  (II),  covering  a  more  limited  range  of  down¬ 
stream  positions,  are  presented  in  Figures  46  to  49  inclusive.  These 
measurements  were  made  at  two  different  values  of  the  jet  outlet  Reynolds 
number  and  velocity  ratio  £. 

An  examination  of  these  velocity  profiles  shows  that  once  the 
ratio  £/Um  exceeds  approximately  .4,  the  tai I  of  the  profi le  constructed 
on  the  basis  of  equation  (2.35)  begins  to  show  a  certain  amount  of  wavi- 
ness,  the  amplitude  of  these  oscillations  or  waves  increasing  as  £/Um 
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Increases.  This  of  course  Is  primarily  due  to  the  fact  that  only  the  first 
four  terms  of  the  asymptotic  expansion  (2.1)  were  used  to  represent  the 
profile.  The  value  of  £/Um  at  which  this  oscillatory  behavior  begins  to 
appear  could  easily  be  Increased  by  simply  Including  more  terms  In  the 
profile  representation.  It  Is  also  Interesting  to  notice  that  the  Inner 
portion  of  the  profile  (i.e.,  T|  <  I)  remains  relatively  well  behaved  for 
all  values  of  £/Um  shown  and  gives  a  very  good  representation  of  the  ex¬ 
perimental  data  in  this  range.  This  is  due  to  the  fact  that  the  similarity 
functions  g  <T|) ,  h(T|),  and  j(7])  as  defined  by  equations  (2.27),  (2.31),  and 
(2.34)  are  extreme  I y  sma II  in  the  region  0  <  T]  <  I .  Thi s  also  bears  out 
the  conclusion  reached  by  Seban  and  Back  (14)  that  the  outer  free  stream 
has  very  little  Influence  on  the  Inner  portion  of  the  flow. 

The  variation  of  6  and  Um  with  the  downstream  coordinate  x  tor 
each  of  the  investigations  mentioned  above  is  presented  In  Figures  50 
through  57.  The  predicted  variations,  as  given  by  equations  (2.48)  and 
(2.49),  are  also  shown.  As  before,  the  constants  C  and  D  were  determined 
to  give  a  reasonable  fit  to  the  experimental  data  in  each  case  and  the 
value  of  ( I4M),  as  defined  by  equation  (2.44),  was  determined  by  simply 
averaging  the  values  of  Um(6F)^^.  Again  the  assumption  that  the  point  of 
zero  shear  stress  coincides  with  the  maximum  velocity  location  (I.e.,  \  = 
1.7977)  was  employed.  It  Is  evident  that  this  assumption  Is  also  reason¬ 
able  in  this  more  general  case  as  far  as  the  prediction  of  the  growth  of 
the  jet  and  decay  of  maximum  velocity  is  concerned.  From  these  figures 
It  is  easily  seen  that  the  variation  of  5/F  with  the  downstream  coordinate 
is  approximately  linear.  This  Is  consistent  with  the  previous  results 
for  the  special  case  £  =  0,  since  in  this  case  the  scale  factor  F  Is  unity. 
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B.  Jet  Outlet  Conditions 

For  a  short  distance  downstream  of  the  outlet,  the  conditions 
in  the  jet  will  be  strongly  dependent  on  the  nature  of  the  flow  immed¬ 
iately  at  the  nozzle.  For  the  plane  turbulent  wall  jet,  the  two  most  im¬ 
portant  factors  which  influence  the  flow  in  this  region,  usually  referred 
to  as  the  jet  core  and  transition  region,  are  the  jet  outlet  Reynolds  num¬ 
ber  and  the  local  geometry  of  the  outlet  itself.  Initial  turbulence  in  the 
jet  has  little  influence  on  the  flow  except  in  a  region  on  the  order  of 
one  slot  width  from  the  jet  outlet. 

Because  of  the  large  amount  of  momentum  in  the  jet,  the  rate 
mass  entrainment  near  the  nozzle  is  very  great  and  the  transverse  velocity 
gradient  i s  at  a  maximum.  This  mixing  process  induces  a  slight  pressure 
rise  in  this  region  and  only  beyond  30  to  40  slot  widths  downstream  of  the 
outlet  does  the  pressure  gradient  become  essentially  zero.  At  this  point 
the  flow  has  already  reached  a  fully  developed  or  self-preserving  state  and 
only  beyond  this  point  is  the  assumption  of  simi  lar  velocity  profi  les  valid. 
Normally,  the  maximum  velocity  in  the  jet  remains  roughly  constant  for  about 
ten  slot  widths  downstream  of  the  outlet  before  decreasing  whi le  the  jet 
thickness  begins  to  increase  immediately.  This  possibly  accounts  for  the 
fact  that  in  many  cases  it  has  been  found  that  the  so-called  virtual  origin 
locations  for  the  maximum  velocity  decay  and  jet  thickness  variation  do 
not  coi nci de. 

In  an  attempt  at  a  more  accurate  determination  of  the  effect  of 
local  conditions  at  the  jet  exit  on  the  downstream  flow,  the  constants  C, 

D,  and  I4W  employed  in  Figures  19  through  26  are  presented  in  Figure  58 
as  functions  of  the  jet  outlet  Reynolds  number.  There  appears  to  be  some 
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evldence  that  each  of  these  factors  is  a  rather  slowly  varying  function 
of  the  Reynolds  number  and  that  the  effect  of  local  geometry  differences 
is  not  overly  great. 

In  addition  to  the  jet  outlet  Reynolds  number  and  local  geome¬ 
trical  configuration,  the  conditions  near  the  jet  outlet  in  the  more  general 
case  of  a  turbulent  wall  jet  beneath  a  secondary  stream  are  also  affected 
by  the  presence  of  the  secondary  stream  itself  and  the  boundary  layer  on 
the  lip  of  the  jet  arising  as  a  result  of  this  secondary  stream.  Much  of 
the  practical  interest  in  wall  jets  is  in  flows  with  an  external  stream 
and  particularly  in  the  mixing  of  a  wall  jet  with  such  a  boundary  layer. 
However,  it  has  been  found  that  except  for  relatively  large  values  of  the 
velocity  ratio  Q,  the  initial  boundary  layer  has  very  little  effect  on 
the  fully  developed  velocity  profiles  which  appear  at  approximately  30 
slot  widths  downstream  of  the  jet  outlet.  That  is,  the  momentum  deficit 
of  the  initial  boundary  layer  is  absorbed  so  quickly  that  the  influence 
of  the  initial  boundary  layer  is  confined  to  the  region  immediately  down¬ 
stream  of  the  outlet.  Specifically,  Bradshaw  and  Gee  (3)  have  stated  that 
the  fully  developed  jet  is  not  greatly  affected  by  the  initial  boundary 
layer  provided  the  jet  momentum  excess  is  approximately  five  times  as  large 
as  the  boundary  layer  momentum  deficit.  Therefore,  for  all  practical  ap¬ 
plications  the  initial  boundary  layer  is  important  only  in  so  far  as  it 
affects  conditions  immediately  downstream  of  the  jet  outlet. 

Qualitatively,  the  presence  of  the  secondary  free  stream  and 
resulting  boundary  layer  above  the  turbulent  wall  jet  affects  the  condi¬ 
tions  in  the  vicinity  of  the  jet  outlet  very  little.  There  is  still  the 
slight  induced  pressure  rise  just  down-stream  of  the  jet  outlet  with  an 
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essential  ly  zero  pressure  gradient  beyond  about  30  to  40  slot  widths  down¬ 
stream  of  the  outlet.  The  potential  core  or  region  of  undiminished  total 
pressure  in  which  the  maximum  velocity  remains  roughly  constant  still  exists 
followed  by  the  transition  region  in  which  the  velocity  profi les  reach  the  - 
fully  developed  state.  Only  beyond  this  point  where  the  velocity  profiles 
become  fully  developed  is  the  expansion  in  terms  of  simi larity  functions 
valid. 

For  this  more  general  case,  the  constants  C,  D,  and  I4M  employed 
in  Figures  50  through  57  are  presented  in  Figure  59  as  functions  of  the 
jet  outlet  Reynolds  number.  It  appears  that  any  dependence  on  the  Reynolds 
number  which  may  exist  has  been  completely  overshadowed  by  the  effects  of 
the  secondary  stream  and,  in  particular,  the  effects  of  the  initial  boun¬ 
dary  layer.  It  is  interesting  to  note  that  the  numerical  range  of  each 
of  these  constants  has  not  been  greatly  affected  by  the  addition  of  the 


secondary  stream. 
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Z  CONCLUSIONS 

1+  has  been  found  that  the  velocity  profi les  in  the  plane  tui — 

bulent  wall  jet  and  turbulent  wall  jet  beneath  a  secondary  uniform  stream 

can  be  adequately  predicted  by  an  asymptotic  expansion  of  the  form 

li  =  fCTp  +  £  g  (Tj)  +£f.  h(T))+£i  j(T)> 

Um  Um  Um3  Um3 

(provided  0  <  (£/Um)  <  .4,)  where 

T]  =  y/6,  C  =  U/Vj 

This  upper  bound  on  the  expansion  parameter  £/Um  could  be  increased  some¬ 
what  by  simply  adding  more  terms  to  the  above  series  expansion. 

Based  on  the  assumption  that  the  shear  stress  distribution  in 
the  jet  can  also  be  represented  by  an  expansion  in  terms  of  similarity 
functions,  the  relation 

Um(6F)l/X  =  constant  =  I  +  M 

between  the  jet  thickness  6  and  maximum  velocity  in  the  profi le  Um  is  ob¬ 
tained.  The  scale  factor  F  is  related  to  the  expansion  parameter  £/Um 
by  the  expression 

1  (  I  +1  )  =  f (F)  +  L  g(F)  +£l  h (F)  +  j£  j (F) 

2  Um  Um  Um3  Um3 

By  considering  the  effect  of  introducing  a  wall  into  a  two- 
dimensional,  free  turbulent  jet,  an  approximate  relation  for  the  growth 
of  the  jet  is  obtained,  namely 

6  -  0  F  exp[C  x(2"X)/X  ] 

where  C  and  0  are  constants.  The  simi larity  parameter  \  was  found  to  be 
dependent  on  the  location  of  the  point  of  zero  shear  stress,  which  lies 
somewhere  in  the  vicinity  of  the  point  of  maximum  velocity  in  the  profi  le. 
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It  was  found  that  the  assumption  that  the  point  of  zero  shear  stress  coin¬ 
cided  with  the  maximum  velocity  location  was  quite  reasonable  as  far  as  the 
prediction  of  the  jet  thickness  and  maximum  velocity  decay  is  concerned. 

Because  of  the  complex  nature  of  the  turbulent  wall  jet  flow, 
exact  solutions  for  the  similarity  functions  f(T|),  g(1p,  h(T|),  and  j(T]) 
could  not  be  obtained.  Instead,  empirical  analytic  expressions  satisfying 
the  required  boundary  conditions  were  adopted  and  these  proved  to  be  quite 
adequate.  Because  the  functions  g(T|),  hCT]),  and  j (T|) ,  as  well  as  any  ad¬ 
ditional  similarity  functions  which  may  be  employed,  are  extremely  small 
in  the  region  0  <  T]  <  I,  the  scale  factor  F  can  be  determined  quite  ac¬ 
curately  by  using  only  very  few  terms  of  the  above  series  expansion.  For 
this  reason,  the  above  relations  for  the  growth  of  the  jet  and  decay  of 
maximum  velocity  are  reasonably  accurate  over  a  much  larger  range  of  the 
expansion  parameter  C/'Um  than  that  for  the  velocity  profiles  themselves. 
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HI  APPENDIX 


I .  Si nau lar i tv  in  f  1  (Tp  : 

I f  the  conditions 

f(l>  =  I 


and 


f'(l)  =  0 

are  used  to  determine  the  constants  A  and  B,  then 

A  =  [I  -  erf (B) ]_1  =  1.14406 
B  *  . I  12060 


and 

fCV1)  =  -5 

Integrating  by  parts,  using  the  new  values  of  A  and  B 

1  (T|)dT)  =  A  Tll/7  [I  -  erf  (BTp  ]dT| 

=  1  +  7  AB  f1  Tl8/7  e~®S^2dTl 
8  4^/ff 

Expanding  in  series  and  integrating  term  by  term 

r1  f  <T|)  dn  =  i  +  49  ab  [  _l  -  ai  +&!-••■] 

0  8  4/rT  1 5  29  86 

=  .933688 

Therefore,  for  the  function  satisfying  conditions  (1.7c)  and  (l.ll). 

Tin 

J  f  (TJ)  dT)  =  .933688  •  T|m  =  .154427 
°o 

T|  T| 

The  numerical  evaluation  of  f  f(z)dz  and  f  f3(z)dz  was  then 

u0  «o 

carried  out  using  a  stepwise  integration  technique  in  which  a  finite  va 

of  f'(0)  was  chosen  such  that 

T| 

Jom  fCTpdTl  =  .154427 
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The  value  of  J*°  f (T|) dT|  was  found  to  be  1.09790  which  can  also  be  determined 
analytically  as  a  further  check  on  the  accuracy  of  the  integration  tech¬ 
nique.  Integrating  by  parts 

J  Alf/7  [I  -  erf (BT]> ]dT|  =  7  AB  j“  lf/7  e-63^2 dT| 

0  4  -/FT 


=  7  AB~8^7  f°  z1^14  e“z  dz 

8/F  00 

where  the  substitution  z  =  BST1®  has  been  made.  This  is  expressable  in 

terms  of  a  Gamma  function  so  that 

f  f  (Tl)dTl  =  7  AB~e/7  T  (15/14)  =  1.09797 
°  4  yn 

which  agrees  almost  exactly  with  the  value  obtained  from  the  stepwise 
i ntegration. 


2.  Example  Velocity  Profile  Construction 

Given  £,  x,  and  the  constants  C,  D,  and  l+M,  the  values  of  6/F 
and  UmF3^  can  be  obtained  from  equations  (2.48)  and  (2.49).  If  a  value 
of  Um  is  chosen,  the  value  of  the  scale  factor  F  as  a  function  of  £/Um 
can  be  determined  from  Figure  4.  A  new  estimate  of  Um  can  then  be  obtained 
since  is  known.  This  procedure  is  then  repeated  until  Um  and  F 

have  been  determined.  The  jet  thickness  6  can  now  be  determined  since 
6/F  i s  known. 

Knowing  £/Um,  values  of  u/Um  can  be  calculated  using  equation 
(2.35)  and  Tables  I,  II,  DU,  and  JJ.  The  non-dimensional  vertical  co¬ 
ordinate  y/5  is  then  obtained  by  scaling  the  values  of  T|  in  the  tables 
by  the  factor  F"1  found  above.  A  plot  of  u/Um  versus  y/6  then  completes 
the  velocity  profi le  construction. 


TABLE  I 


V 

«1>) 

/fdt, 

S(i>) 

3'(7J) 

.000 

.00000 

- 

.0000 

.00000 

1.00000 

.0000 

.010 

.74040 

10.2770 

•  0047 

.00460 

.96738 

-4.1466 

.020 

.89308 

9.3039 

.0147 

.01101 

.92147 

-4.9338 

.030 

.87991 

3.4899 

.0232 

.01834 

.86929 

-3.4497 

.040 

.90392 

2.9204 

.0322 

.02629 

.81297 

-3.7933 

.090 

.92748 

1.9134 

.0413 

.03470 

.73368 

-6.0489 

.040 

.94440 

1.4931 

.0307 

.04347 

.69220 

-6.2387 

.070 

.93771 

1.1827 

.0402 

.09232 

.62906 

-6.3811 

.080 

.94829 

.9430 

.0498 

.06180 

.36470 

-6.4866 

.090 

.97473 

.7313 

.0794 

.07126 

•49943 

-6.9624 

.100 

.98343 

.3944 

.0894 

.08086 

.43333 

-6.6138 

.120 

.99277 

.3919 

.1091 

.10041 

•30070 

-6.6577 

.140 

.99792 

.1719 

.1290 

.12023 

.16763 

-6.6401 

.140 

.99991 

.0324 

•  1490 

.14020 

.03941 

-6.3735 

.180 

.99941 

-.0790 

.1690 

.16019 

-.09314 

-6.4758 

.200 

.99488 

-.1703 

.1890 

.18013 

-.22334 

-6.5424 

.300 

.94383 

-.4948 

.2873 

.27670 

-.81368 

-5.4033 

.400 

.91030 

-.9948 

.3811 

.36478 

-1.29683 

-4.2425 

.300 

.84447 

-.4714 

.4690 

.44211 

-1.66124 

-3.0507 

.400 

.77793 

-.7041 

.3903 

.30813 

-1.90940 

-1.9301 

.700 

.70433 

-.7142 

.6249 

.96324 

-2.09137 

-.9379 

.800 

.43332 

-.7032 

.6919 

.60827 

-2.10224 

-.1031 

.900 

.94419 

-.4783 

.7916 

•64438 

-2.07779 

.3642 

1.000 

.30000 

-.4432 

.8049 

.67280 

-1.99487 

1.0677 

1.100 

.43774 

-.4010 

.8917 

.69478 

-1.86931 

1.4198 

1.200 

.37999 

-.9341 

.<926 

.71149 

-1.71333 

1.6389 

1.300 

.32700 

-.9044 

.9279 

.72398 

-1.34329 

1.7462 

1.400 

.27904 

-.4941 

.9981 

.73316 

-1.36910 

1.7636 

1.900 

.23419 

-.4041 

.9838 

•73979 

-1.19480 

1.7123 

1.400 

.19817 

-.3338 

1.0099 

•74490 

-1.02823 

1.6122 

1.700 

.14491 

-.3099 

1.0236 

.74779 

-.87344 

1.4796 

1.800 

.13409 

-.2472 

1.0387 

.79003 

-.73292 

1.3288 

1.900 

.11134 

-.2280 

1.0910 

.73198 

-.60791 

1.1710 

2.000 

.09034 

-.1924 

1.0610 

.79299 

-.49867 

1.0146 

2.100 

.07270 

-.1411 

1.0692 

.79326 

-.40474 

.8637 

2.200 

.03801 

-.1333 

1.0797 

.79368 

-.32316 

.7282 

2.300 

.04989 

-.1099 

1.0809 

.73393 

-.29864 

.6046 

2.400 

.03999 

-.0889 

1.0849 

.79412 

-.20374 

.4938 

2.900 

.02799 

-.0714 

1.0881 

.73422 

-.13898 

.4018 

2.400 

.02199 

-.0970 

1.0906 

.73428 

-.12291 

.3220 

2.700 

.01430 

-.0430 

1.0929 

.73432 

-.09419 

.2992 

2.800 

.01291 

-.0392 

1.0939 

.79434 

-.07147 

.2002 

2.900 

.00940 

-.0273 

1.0990 

.79433 

-.09376 

.1594 

9.000 

.00700 

-.0209 

1.0998 

.79439 

-.04009 

.1199 

TABLE  ft 


fltol 

_ 1*  adw _ 

_ rial _ 

_ Hal _ 

- - —  —  *b 

<000 

.00000 

.0000 

.0000 

.0000 

.000 

.010 

.00000 

.0000 

.0000 

.1785 

22.691 

.020 

.00000 

.0000 

•oooo 

.9297 

27.108 

.0)0 

.00000 

.0000 

.0000 

.7133 

29.621 

.0*0 

.00000 

.0000 

.0000 

1.0233 

31.71* 

.0)0 

.00000 

.0000 

.0000 

1.3979 

33.101 

.060 

.00000 

.0000 

.0000 

1.6899 

39.139 

.070 

.00000 

.0000 

.0000 

2.0298 

39.918 

.080 

.00000 

.0000 

.0000 

2.3821 

35.996 

.090 

•00000 

.0000 

.0000 

2.7392 

35.911 

.100 

.00000 

.0000 

.0000 

3.0998 

36.192 

.120 

.00000 

.0000 

.0000 

3.8267 

36.932 

.190 

.00000 

.0000 

.0000 

9.3399 

36.336 

.160 

.00000 

.0000 

•  0000 

5.278* 

35.982 

.180 

.00000 

.0000 

.0000 

3.9928 

35.926 

.200 

.00000 

.0000 

.0000 

6.6999 

39.707 

.300 

.00000 

.0000 

.0000 

9.9298 

29.568 

.900 

.00000 

.0000 

.0000 

12.3687 

23.216 

.300 

.00000 

.0000 

.0000 

19.9628 

16.699 

.600 

.00000 

.0000 

.0000 

13.9208 

10.562 

.700 

.00000 

.0000 

.0000 

16.6988 

3.132 

.800 

.00000 

.0000 

.0000 

16.9760 

.56* 

.900 

.00001 

.0002 

.0000 

16.8922 

-3.088 

1.000 

.00006 

.0010 

.0000 

16.3883 

-3.8*7 

1.100 

.00029 

.0090 

•  0000 

13.7001 

-7.788 

1.200 

.00109 

•0129 

.0001 

19.8590 

-9.029 

1.300 

.00313 

.0322 

•  0003 

13.9128 

-9.717 

1.900 

.00811 

.0709 

•  0008 

12.9236 

-10.017 

1.500 

.01821 

.1360 

.0021 

11.9171 

-10.089 

1.600 

.03631 

.2308 

.0097 

10.9091 

-10.063 

1.700 

.06)26 

.3317 

.0097 

9.9053 

-10.013 

1.800 

.10716 

.987* 

.0182 

0*9066 

-9.935 

1.900 

.16268 

.6212 

.0316 

7.9136 

-9.83? 

2.000 

.23072 

.7330 

.0512 

6.9394 

-9.6*8 

2.100 

.30839 

.8197 

.0780 

3.9912 

-9.288 

2.200 

.39223 

.832) 

.1131 

3.0879 

-8.798 

2.300 

.97738 

.897) 

.1363 

9.2979 

-8.037 

2.900 

.36032 

.8037 

.2083 

3.9899 

-7.193 

2.900 

.63781 

.7369 

.2689 

2.8109 

-6.277 

2.600 

.70729 

.6301 

.3358 

2.2301 

-3.392 

2.700 

.76760 

.3367 

.9096 

1.7919 

-9.9*1 

2.800 

.81861 

.96*0 

.9890 

1.3393 

-3.613 

2.900 

.86062 

.3777 

.3730 

1.0137 

-2.881 

3.000 

.89996 

.3009 

.6608 

.7398 

-2.233 

TABLE  UL 


h(«) ..  - 

h'fwl 

_ Lhdfi _ 

WLai _ 

-r  —  %  ' 

.000 

.00000 

.0000 

.0000 

.000 

.00 

.010 

.00000 

.0000 

.0000 

-.939 

-119.29 

.020 

.00000 

.0000 

.0000 

-2.239 

-142.32 

.030 

.00000 

.0000 

.0000 

-3.761 

-136,79 

.040 

.00000 

.0000 

.0000 

-3.381 

-166.73 

.030 

•00000 

.0000 

.0000 

-7.086 

-174.02 

.060 

.00000 

.0000 

•  0000 

-8.833 

-179.48 

.070 

*00000 

•  0000 

.0000 

-10.672 

-103*50 

.000 

.00000 

.0000 

•  0000 

-12.523 

-186.62 

.090 

•00000 

.0000 

.0000 

-14.401 

-188.80 

.100 

.00000 

.0000 

.0000 

-16.297 

-190.28 

.120 

•00000 

.0000 

•  0000 

-20.118 

-191.54 

.140 

.00000 

•  0000 

•  0000 

-23.947 

-191.03 

.160 

•00000 

•  0000 

.0000 

-27.731 

-189.17 

.180 

.00000 

•  0000 

•  0000 

-31*506 

-186.25 

.200 

.00000 

•  0000 

.0000 

-35.193 

-182.47 

.300 

•00000 

•  0000 

.0000 

-52.178 

-155.45 

.400 

•00000 

.0000 

.0000 

-66.070 

-122.06 

.300 

.00000 

.0000 

•  oooo 

-76.562 

-87.77 

.600 

.ooooo 

.0000 

.0000 

-83.701 

-53.53 

.700 

.00000 

.0001 

.0000 

-87.792 

-26.98 

.800 

.00004 

.0009 

.0000 

-89.250 

-2.07 

•  900 

.00030 

.0034 

.0000 

-88.347 

16.22 

1.000 

.00137 

.0237 

.0001 

-86.165 

30*64 

1.100 

.00623 

.0787 

.0004 

-82.568 

40.60 

1.200 

.01977 

.2073 

.0016 

-78.180 

46.53 

1.300 

.03144 

•  4466 

.0050 

-73.377 

49.00 

1*400 

.11313 

,8035 

.0129 

-68.469 

48.81 

1.300 

.21507 

1.2380 

.0290 

-63.668 

47.04 

1*600 

.33980 

1.6411 

.0574 

-59.072 

44.91 

1.700 

.33808 

1.8890 

.1021 

-54.662 

43.44 

1*000 

.72912 

1.8852 

.1654 

-30.329 

43.39 

1.900 

.90579 

1.6041 

.2474 

-45.942 

44*54 

2.000 

1.04241 

1.0988 

.3432 

-41.402 

46.29 

2.100 

1.12161 

.4759 

.4340 

-36.694 

47.74 

2.200 

1.13765 

-.1437 

.3674 

-31.891 

48.08 

2.300 

1.09577 

-.6700 

.6796 

-27.130 

46*06 

2.400 

1*00077 

-1.0427 

.7851 

-22.375 

o 

o 

• 

* 

2.300 

.89274 

-1.2313 

.8803 

-18.374 

39.82 

2.600 

.76333 

-1.3130 

.9632 

-14.639 

34.78 

2.700 

.63340 

-1.2684 

1.0330 

-11.429 

29.41 

2.800 

.51201 

-1.1305 

1 . 0902 

-8.734 

24.14 

2.900 

.40452 

-.9936 

1.1359 

-6.387 

19.30 

3.000 

.31323 

-.8301 

1.1716 

-4,874 

13.06 

TABLE  BE 


_ ILbI _ 

_ ft&l _ 

_ Lida _ 

_ 7(m\ 

_ _ 

- - ~7 - *0  T 

,000 

.00000 

.0000 

.0000 

.000 

.00 

.010 

.00000 

.0000 

.0000 

5.304 

674.20 

•  020 

.00000 

.0000 

.0000 

12.766 

605.45 

.0)0 

.00000 

.0000 

.0000 

21.233 

886.08 

•  040 

.00000 

.0000 

.0000 

30.410 

942.29 

.0)0 

.00000 

.0000 

.0000 
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FIGURE  5.  SECTIONAL  VIEW  OF  WALL  JET 


FIGURE  6.  TOTAL  HEAO  PROBE  AND  CARRIAGE 
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FIGURE  15.  EXPERIMENTAL  JET  THICKNESS  DATA 
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